Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Validation of a metabolite-GWAS network for Populus trichocarpa family 1 UDP-glycosyltransferases.

Metabolite genome-wide association studies (mGWASs) are increasingly used to discover the genetic basis of target phenotypes in plants such as Populus trichocarpa, a biofuel feedstock and model woody plant species. Despite their growing importance in plant genetics and metabolomics, few mGWASs are experimentally validated. Here, we present a functional genomics workflow for validating mGWAS-predicted enzyme-substrate relationships. We focus on uridine diphosphate-glycosyltransferases (UGTs), a large family of enzymes that catalyze sugar transfer to a variety of plant secondary metabolites involved in defense, signaling, and lignification. Glycosylation influences physiological roles, localization within cells and tissues, and metabolic fates of these metabolites. UGTs have substantially expanded in P. trichocarpa, presenting a challenge for large-scale characterization. Using a high-throughput assay, we produced substrate acceptance profiles for 40 previously uncharacterized candidate enzymes. Assays confirmed 10 of 13 leaf mGWAS associations, and a focused metabolite screen demonstrated varying levels of substrate specificity among UGTs. A substrate binding model case study of UGT-23 rationalized observed enzyme activities and mGWAS associations, including glycosylation of trichocarpinene to produce trichocarpin, a major higher-order salicylate in P. trichocarpa. We identified UGTs putatively involved in lignan, flavonoid, salicylate, and phytohormone metabolism, with potential implications for cell wall biosynthesis, nitrogen uptake, and biotic and abiotic stress response that determine sustainable biomass crop production. Our results provide new support for in silico analyses and evidence-based guidance for in vivo functional characterization.

Plant myo-inositol transport influences bacterial colonization phenotypes.

Plant microbiomes are assembled and modified through a complex milieu of biotic and abiotic factors. Despite dynamic and fluctuating contributing variables, specific host metabolites are consistently identified as important mediators of microbial interactions. We combine information from a large-scale metatranscriptomic dataset from natural poplar trees and experimental genetic manipulation assays in seedlings of the model plant Arabidopsis thaliana to converge on a conserved role for transport of the plant metabolite myo-inositol in mediating host-microbe interactions. While microbial catabolism of this compound has been linked to increased host colonization, we identify bacterial phenotypes that occur in both catabolism-dependent and -independent manners, suggesting that myo-inositol may additionally serve as a eukaryotic-derived signaling molecule to modulate microbial activities. Our data suggest host control of this compound and resulting microbial behavior are important mechanisms at play surrounding the host metabolite myo-inositol.

Evaluating the performance of random forest and iterative random forest based methods when applied to gene expression data.

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

Antiviral Strategies Against SARS-CoV-2: A Systems Biology Approach.

The unprecedented scientific achievements in combating the COVID-19 pandemic reflect a global response informed by unprecedented access to data. We now have the ability to rapidly generate a diversity of information on an emerging pathogen and, by using high-performance computing and a systems biology approach, we can mine this wealth of information to understand the complexities of viral pathogenesis and contagion like never before. These efforts will aid in the development of vaccines, antiviral medications, and inform policymakers and clinicians. Here we detail computational protocols developed as SARS-CoV-2 began to spread across the globe. They include pathogen detection, comparative structural proteomics, evolutionary adaptation analysis via network and artificial intelligence methodologies, and multiomic integration. These protocols constitute a core framework on which to build a systems-level infrastructure that can be quickly brought to bear on future pathogens before they evolve into pandemic proportions.

Response to comment on 'SARS-CoV-2 suppresses anticoagulant and fibrinolytic gene expression in the lung'.

Early in the SARS-CoV-2 pandemic, we compared transcriptome data from hospitalized COVID-19 patients and control patients without COVID-19. We found changes in procoagulant and fibrinolytic gene expression in the lungs of COVID-19 patients (Mast et al., 2021). These findings have been challenged based on issues with the samples (Fitzgerald and Jamieson, 2022). We have revisited our previous analyses in the light of this challenge and find that these new analyses support our original conclusions.

Polypharmacy-associated risk of hospitalisation among people ageing with and without HIV: an observational study.

BACKGROUND: Polypharmacy, defined as use of five or more medications concurrently, is associated with adverse health outcomes and people ageing with HIV might be at greater risk than similar uninfected individuals. We aimed to determine whether known pairwise drug interactions (KPDIs) were associated with risk of admission to hospital (hereafter referred to as hospitalisation) and medication count among people ageing with and without HIV after accounting for physiological frailty. METHODS: In this observational study, we collected individual-level data for participants of the Veterans Aging Cohort Study (VACS) with HIV on antiretroviral therapy (ART) and with supressed HIV-1 RNA and people without HIV who were receiving at least one prescription medication, based on active medications in the 2009 fiscal year (ie, Oct 1, 2008, to Sept 30, 2009). We identified KPDIs among these patients by linking prescription fill and refill data with data from DrugBank (version 5.0.11). We collected data on all-cause mortality and hospitalisations between Oct 1, 2009, and March 31, 2019. We compared KPDI counts using random selection and actual patterns of use across medication counts from two to 12. We created a weighted KPDI Index on the basis of the average association of each KPDI with mortality among people ageing without HIV and used nested Cox models stratified by HIV status to estimate the association between medication count and hospitalisation, with incremental adjustments for demographics, physiological frailty, and KPDI Index. FINDINGS: We collected data for 9186 people ageing with HIV and 37 930 individuals without HIV. 45 913 (97·4%) of 47 116 patients were men and the sample was predominantly aged 50-64 years (30 413 [64·6%]). Compared with a random sample of medications, real-world pattern of medication counts and combinations were associated with five-to-six times more KPDIs (eg, for a combination of six medications, KPDI count was 1·09 in the random sample, 5·49 in the HIV-negative population, and 7·13 in the HIV-positive population). For each additional observed medication, people ageing with HIV had approximately 2·94 additional KPDIs and comparators had approximately 2·67 additional KPDIs. Adjustment for demographics, physiological frailty, and KPDI Index reduced the association between medication count and risk of hospitalisation for people ageing with HIV (hazard ratio 1·08 [95% CI 1·07-1·09] reduced to 1·06 [1·05-1·07]) and those without HIV (1·08 [1·07-1·08] reduced to 1·04 [1·03-1·05]). INTERPRETATION: For each additional medication, people ageing with HIV have more drug-drug interactions than those without HIV. Adjusting for known non-ART drug-drug interactions, each additional non-ART medication confers excess risk of hospitalisation for people ageing with HIV. Randomised trials will be needed to determine whether reducing these interactions improves outcomes. FUNDING: National Institutes of Health, National Institute on Alcohol Abuse and Alcoholism, Department of Veterans Affairs Health Services Research & Development, and Office of Research and Development.

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes.

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

Image-Based Methods to Score Fungal Pathogen Symptom Progression and Severity in Excised Arabidopsis Leaves.

Image-based symptom scoring of plant diseases is a powerful tool for associating disease resistance with plant genotypes. Advancements in technology have enabled new imaging and image processing strategies for statistical analysis of time-course experiments. There are several tools available for analyzing symptoms on leaves and fruits of crop plants, but only a few are available for the model plant Arabidopsis thaliana (Arabidopsis). Arabidopsis and the model fungus Botrytis cinerea (Botrytis) comprise a potent model pathosystem for the identification of signaling pathways conferring immunity against this broad host-range necrotrophic fungus. Here, we present two strategies to assess severity and symptom progression of Botrytis infection over time in Arabidopsis leaves. Thus, a pixel classification strategy using color hue values from red-green-blue (RGB) images and a random forest algorithm was used to establish necrotic, chlorotic, and healthy leaf areas. Secondly, using chlorophyll fluorescence (ChlFl) imaging, the maximum quantum yield of photosystem II (Fv/Fm) was determined to define diseased areas and their proportion per total leaf area. Both RGB and ChlFl imaging strategies were employed to track disease progression over time. This has provided a robust and sensitive method for detecting sensitive or resistant genetic backgrounds. A full methodological workflow, from plant culture to data analysis, is described.

Potential Pathogenicity Determinants Identified from Structural Proteomics of SARS-CoV and SARS-CoV-2.

Despite SARS-CoV and SARS-CoV-2 being equipped with highly similar protein arsenals, the corresponding zoonoses have spread among humans at extremely different rates. The specific characteristics of these viruses that led to such distinct outcomes remain unclear. Here, we apply proteome-wide comparative structural analysis aiming to identify the unique molecular elements in the SARS-CoV-2 proteome that may explain the differing consequences. By combining protein modeling and molecular dynamics simulations, we suggest nonconservative substitutions in functional regions of the spike glycoprotein (S), nsp1, and nsp3 that are contributing to differences in virulence. Particularly, we explain why the substitutions at the receptor-binding domain of S affect the structure-dynamics behavior in complexes with putative host receptors. Conservation of functional protein regions within the two taxa is also noteworthy. We suggest that the highly conserved main protease, nsp5, of SARS-CoV and SARS-CoV-2 is part of their mechanism of circumventing the host interferon antiviral response. Overall, most substitutions occur on the protein surfaces and may be modulating their antigenic properties and interactions with other macromolecules. Our results imply that the striking difference in the pervasiveness of SARS-CoV-2 and SARS-CoV among humans seems to significantly derive from molecular features that modulate the efficiency of viral particles in entering the host cells and blocking the host immune response.

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models.

BACKGROUND: A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic. RESULTS: Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus. CONCLUSIONS: These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa.

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.

Network Modeling of Complex Data Sets.

We demonstrate a selection of network and machine learning techniques useful in the analysis of complex datasets, including 2-way similarity networks, Markov clustering, enrichment statistical networks, FCROS differential analysis, and random forests. We demonstrate each of these techniques on the Populus trichocarpa gene expression atlas.

Neurotransmitter networks in mouse prefrontal cortex are reconfigured by isoflurane anesthesia.

This study quantified eight small-molecule neurotransmitters collected simultaneously from prefrontal cortex of C57BL/6J mice (n = 23) during wakefulness and during isoflurane anesthesia (1.3%). Using isoflurane anesthesia as an independent variable enabled evaluation of the hypothesis that isoflurane anesthesia differentially alters concentrations of multiple neurotransmitters and their interactions. Machine learning was applied to reveal higher order interactions among neurotransmitters. Using a between-subjects design, microdialysis was performed during wakefulness and during anesthesia. Concentrations (nM) of acetylcholine, adenosine, dopamine, GABA, glutamate, histamine, norepinephrine, and serotonin in the dialysis samples are reported (means ± SD). Relative to wakefulness, acetylcholine concentration was lower during isoflurane anesthesia (1.254 ± 1.118 vs. 0.401 ± 0.134, P = 0.009), and concentrations of adenosine (29.456 ± 29.756 vs. 101.321 ± 38.603, P < 0.001), dopamine (0.0578 ± 0.0384 vs. 0.113 ± 0.084, P = 0.036), and norepinephrine (0.126 ± 0.080 vs. 0.219 ± 0.066, P = 0.010) were higher during anesthesia. Isoflurane reconfigured neurotransmitter interactions in prefrontal cortex, and the state of isoflurane anesthesia was reliably predicted by prefrontal cortex concentrations of adenosine, norepinephrine, and acetylcholine. A novel finding to emerge from machine learning analyses is that neurotransmitter concentration profiles in mouse prefrontal cortex undergo functional reconfiguration during isoflurane anesthesia. Adenosine, norepinephrine, and acetylcholine showed high feature importance, supporting the interpretation that interactions among these three transmitters may play a key role in modulating levels of cortical and behavioral arousal.NEW & NOTEWORTHY This study discovered that interactions between neurotransmitters in mouse prefrontal cortex were altered during isoflurane anesthesia relative to wakefulness. Machine learning further demonstrated that, relative to wakefulness, higher order interactions among neurotransmitters were disrupted during isoflurane administration. These findings extend to the neurochemical domain the concept that anesthetic-induced loss of wakefulness results from a disruption of neural network connectivity.

Plasticity of Phymatotrichopsis omnivora infection strategies is dependent on host and nonhost plant responses.

Necrotrophic fungi constitute the largest group of plant fungal pathogens that cause heavy crop losses worldwide. Phymatotrichopsis omnivora is a broad host, soil-borne necrotrophic fungal pathogen that infects over 2,000 dicotyledonous plants. The molecular basis of such broad host range is unknown. We conducted cell biology and transcriptomic studies in Medicago truncatula (susceptible), Brachypodium distachyon (resistant/nonhost), and Arabidopsis thaliana (partially resistant) to understand P. omnivora virulence mechanisms. We performed defence gene analysis, gene enrichments, and correlational network studies during key infection stages. We identified that P. omnivora infects the susceptible plant as a traditional necrotroph. However, it infects the partially resistant plant as a hemi-biotroph triggering salicylic acid-mediated defence pathways in the plant. Further, the infection strategy in partially resistant plants is determined by the host responses during early infection stages. Mutant analyses in A. thaliana established the role of small peptides PEP1 and PEP2 in defence against P. omnivora. The resistant/nonhost B. distachyon triggered stress responses involving sugars and aromatic acids. Bdwat1 mutant analysis identified the role of cell walls in defence. This is the first report that describes the plasticity in infection strategies of P. omnivora providing insights into broad host range.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence.

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

Plant Host-Associated Mechanisms for Microbial Selection.

Plants serve as host to numerous microorganisms. The members of these microbial communities interact among each other and with the plant, and there is increasing evidence to suggest that the microbial community may promote plant growth, improve drought tolerance, facilitate pathogen defense and even assist in environmental remediation. Therefore, it is important to better understand the mechanisms that influence the composition and structure of microbial communities, and what role the host may play in the recruitment and control of its microbiome. In particular, there is a growing body of research to suggest that plant defense systems not only provide a layer of protection against pathogens but may also actively manage the composition of the overall microbiome. In this review, we provide an overview of the current research into mechanisms employed by the plant host to select for and control its microbiome. We specifically review recent research that expands upon the role of keystone microbial species, phytohormones, and abiotic stress, and in how they relate to plant driven dynamic microbial structuring.

Multi-Phenotype Association Decomposition: Unraveling Complex Gene-Phenotype Relationships.

Various patterns of multi-phenotype associations (MPAs) exist in the results of Genome Wide Association Studies (GWAS) involving different topologies of single nucleotide polymorphism (SNP)-phenotype associations. These can provide interesting information about the different impacts of a gene on closely related phenotypes or disparate phenotypes (pleiotropy). In this work we present MPA Decomposition, a new network-based approach which decomposes the results of a multi-phenotype GWAS study into three bipartite networks, which, when used together, unravel the multi-phenotype signatures of genes on a genome-wide scale. The decomposition involves the construction of a phenotype powerset space, and subsequent mapping of genes into this new space. Clustering of genes in this powerset space groups genes based on their detailed MPA signatures. We show that this method allows us to find multiple different MPA and pleiotropic signatures within individual genes and to classify and cluster genes based on these SNP-phenotype association topologies. We demonstrate the use of this approach on a GWAS analysis of a large population of 882 Populus trichocarpa genotypes using untargeted metabolomics phenotypes. This method should prove invaluable in the interpretation of large GWAS datasets and aid in future synthetic biology efforts designed to optimize phenotypes of interest.

Fungal-Bacterial Networks in the Populus Rhizobiome Are Impacted by Soil Properties and Host Genotype.

Plant root-associated microbial symbionts comprise the plant rhizobiome. These microbes function in provisioning nutrients and water to their hosts, impacting plant health and disease. The plant microbiome is shaped by plant species, plant genotype, soil and environmental conditions, but the contributions of these variables are hard to disentangle from each other in natural systems. We used bioassay common garden experiments to decouple plant genotype and soil property impacts on fungal and bacterial community structure in the Populus rhizobiome. High throughput amplification and sequencing of 16S, ITS, 28S and 18S rDNA was accomplished through 454 pyrosequencing. Co-association patterns of fungal and bacterial taxa were assessed with 16S and ITS datasets. Community bipartite fungal-bacterial networks and PERMANOVA results attribute significant difference in fungal or bacterial communities to soil origin, soil chemical properties and plant genotype. Indicator species analysis identified a common set of root bacteria as well as endophytic and ectomycorrhizal fungi associated with Populus in different soils. However, no single taxon, or consortium of microbes, was indicative of a particular Populus genotype. Fungal-bacterial networks were over-represented in arbuscular mycorrhizal, endophytic, and ectomycorrhizal fungi, as well as bacteria belonging to the orders Rhizobiales, Chitinophagales, Cytophagales, and Burkholderiales. These results demonstrate the importance of soil and plant genotype on fungal-bacterial networks in the belowground plant microbiome.

Epigenetic effects of paternal cocaine on reward stimulus behavior and accumbens gene expression in mice.

Paternal cocaine use causes phenotypic alterations in offspring behavior and associated neural processing. In rodents, changes in first generation (F1) offspring include drug reward behavior, circadian timing, and anxiety responses. This study, utilizing a murine (C57BL/6J) oral cocaine model, examines the effects of paternal cocaine exposure on fundamental characteristics of offspring reward responses, including: 1) the extent of cocaine-induced effects after different durations of sire drug withdrawal; 2) sex- and drug-dependent differences in F1 reward preference; 3) effects on second generation (F2) cocaine preference; and 4) corresponding changes in reward area (nucleus accumbens) mRNA expression. We demonstrate that paternal cocaine intake over a single ˜40-day spermatogenic cycle significantly decreased cocaine (but not ethanol or sucrose) preference in a sex-specific manner in F1 mice from sires mated 24 h after drug withdrawal. However, F1 offspring of sires bred 4 months after withdrawal did not exhibit altered cocaine preference. Altered cocaine preference also was not observed in F2's. RNASeq analyses of F1 accumbens tissue revealed changes in gene expression in male offspring of cocaine-exposed sires, including many genes not previously linked to cocaine addiction. Enrichment analyses highlight genes linked to CNS development, synaptic signaling, extracellular matrix, and immune function. Expression correlation analyses identified a novel target, Fam19a4, that may negatively regulate many genes in the accumbens, including genes already identified in addiction. Collectively, these results reveal that paternal cocaine effects in F1 offspring may involve temporally limited epigenetic germline effects and identify new genetic targets for addiction research.